Astronomy & Astrophysics manuscript no. arxiv'version'polarisation'201 1 


©ESO 2011 


October 25, 2011 





Spatially resolved hard X-ray polarization in solar flares: effects of 
Compton scattering and bremsstrahlung 



N. L. S. Jeffrey and E. P. Kontar 



o 

(N 
-i— > 
O 

O 

(N 
& 

6 



> 

m 

OS 

o> 
d 



X 



School of Physics & Astronomy, University of Glasgow, G12 8QQ, United Kingdom 
Received ; Accepted 

ABSTRACT 

Aims. To study the polarization of hard X-ray (HXR) sources in the solar atmosphere, including Compton backscattering of photons 
in the photosphere (the albedo effect) and the spatial distribution of polarization across the source. 

Methods. HXR photon polarization and spectra produced via electron-ion bremsstrahlung emission are calculated from various elec- 
tron distributions typical for solar flares. Compton scattering and photoelectric absorption are then modelled using Monte Carlo 
simulations of photon transport in the photosphere to study the observed (primary and albedo) sources. Polarization maps across 
HXR sources (primary and albedo components) for each of the modelled electron distributions are calculated at various source loca- 
tions from the solar centre to the limb. 

Results. We show that Compton scattering produces a distinct polarization variation across the albedo patch at peak albedo energies 
of 20-50 keV for all anisotropies modelled. The results show that there are distinct spatial polarization changes in both the radial and 
perpendicular to radial directions across the extent of the HXR source at a given disk location. In the radial direction, the polarization 
magnitude and direction at specific positions along the HXR source will either increase or decrease with increased photon distribution 
directivity towards the photosphere. We also show how high electron cutoff energies influence the direction of polarization at above 
~100keV. 

Conclusions. Spatially resolved HXR polarization measurements can provide important information about the directivity and en- 
ergetics of the electron distribution. Our results indicate the preferred angular resolution of polarization measurements required to 
distinguish between the scattered and primary components. We also show how spatially resolved polarization measurements could be 
used to probe the emission pattern of an HXR source, using both the magnitude and the direction of the polarization. 

Key words. Sun: Flares - Sun: X-rays, gamma rays - Sun: corona - Sun: chromosphere - Sun: polarization/polarisation 



1. Introduction 

Solar flare accelerated electrons can be studied via their radia- 
tion using various emission mechanisms. Electrons with ener- 
gies from tens to hundreds of keV can reach the denser lay- 
ers of the solar chromosphere and produce strong hard X-ray 
(HXR) emission via bremsstrahlung; the properties of which 
tell us about the energetics and directivity of the accelerated 
electron distribution. A major insight into the angular prop- 
erties of the photon distribution and hence the radiating elec- 
tron distribution can come directly from the polarization of the 
HXRs. Since the anisotropy and polarization of the photon dis- 
tribution produced by bremsstrahlung wil l increase with the 
aniso t ropy of the electron dis tribution (e.g. lBrown|[T972t lHaugl 
1 19721 iLeach & Petrosianlll983h . then in theory HXR polariza- 
tion is one of the direct tools available to infer the anisotropy 
of the emitting electron distribution. However, compared with 
the other HXR observables: energy, spatial location, source size, 
time of ar rival etc., pol arization measurements through the years 
(see lKontar et al.l201 ll as a review) have been fraught with diffi- 
culties and the measurements obtained have often been met with 
skepticism. Nonetheless many missions have reported measure- 
ments of HXR polarization from solar flares (Tindo et al. 197( 
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19721 iNakada et al]l 19741 iTindo et al J[l976k iLemen et alJl!982 
McConne ll et alJ l2004t iBoggs et all 120061: ISuarez-Garcia et alJ 
2006tlMcConnell et al]l2007l) . 
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Another valuable tool for determining the directivity of the 
photon distribution is the Compton backscattered (albedo) com- 
ponent. Photons emitted towards the photosphere can Compton 
scatter back out towards the observer and be detected alongside 
the primary photons directly emitted from the source (Tomblin 
[19721 ISantaneelo et al J 19731) . Even for an isotropic source (min- 
imum albedo), the albedo component can account for as much as 
40% of the observed flux at peak albedo energies of 20-50 keV 
dBai & Ramatvll978tlKontar et al.l2006l) and for high downward 
beaming, albedo can dominate the observed flux at certain ener- 
gies. Separating the albedo component from the directly emitted 
primary component is essential as it provides us with a unique 
(and only for individual flares) opportunity of inferring the prop- 
erties of the photon distribution directed towards th e photosphere 
dKontar & Brownll2006l: iBattaglia & Kontarll201 lb . which other- 
wise would be poorly known. The albedo component can also 
change various measurable properties of HXR sources, from 
producing a dist inctive change in the photon spectrum between 
~ 10 - 100 keV dBai & Ramatvll 19781) to increasing source size s 
and shifting the position of sources dKontar &.Teffrevl l2010). 
Compton scattering is also polarization dependent and hence the 
polariz ation of an HXR source will be altered by the albedo com- 
ponen t dHenouxl 19751: lLanger & Petrosianlll977t iBai & Ra matv 
1978). The albedo component will not only change the polariza- 
tion of the entire source but it will also produce spatial variations 
in polarization across the extent of the albedo source. 

Both HXR polarization and albedo measurements contribute 
information regarding the directivity of the photon distribution 
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and hence allow us to deduce the pitc h angle dis tribution of the 
X-ray emitting electrons (e.g. iKontar et alJl201 lL as a review). 
In addition, the direction of the polarization vector is related to 
the direction of electron bea ming and hence the dominant di- 
rection of electrons in a loop (lEmslie et alj|2008h . Spatially re- 
solved polarization measurements across an HXR source caused 
by albedo, have the advantage over spatially integrated measure- 
ments because both the magnitude and the direction of polariza- 
tion will change with photon directivity and allow us to map the 
albedo and primary components. Therefore understanding how 
these two parameters change with the photon anisotropy is es- 
sential and provides us with a new method of investigating the 
entire photon anisotropy from a single HXR source. 

Upcoming missions such as the G amma-Ray 
Imag er/Polarimeter for Solar flares (GRIPS) ( Shi h et all 
2009) (which is planned to be flown on a balloon in 2012) or 
other astrophysical missions (Bloseret ail I2009L iMuleri et al. 
2009) may be able to provide us with spatially resolved po- 
larization measurements. Though GRIPS will have an angular 
resolution of ~ 12" available at energies as low as ~ 25 
keV, it is doubtful that the GRIPS polarimeter will be able to 
measure polarization at peak albedo energies of 20-50 keV, 
due to photoelectric absorption dominating over Compton 
scattering in the detector at these energies, but hopefully some 
technologies involved in GRIPS can be used for future missions. 
Nevertheless, although there are a number of sim ulations for the 
spatia l ly integrated polarization signal in flares (lElwert & Haug 
19701: iHaud Il972h lLanger & Petrosianl 1 1 977h iB ai & Ramatvl 
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19781 IZharkova et alJll995t lEmslie etallbooa IZharkova et alJ 
2010), there are currently no theoretical predictions or modelled 
spatially resolved hard X-ray polarization maps. In addition, 
the effect of X-ray albedo on polarization signatures is often 
ignored in the models. 

In this paper, we compute spatially resolved polarization 
across HXR sources at various locations on the solar disk, taking 
into account the influence of albedo, for various emitting elec- 
tron populations. We also examine how changing the maximum 
electron energy available during bremsstrahlung can alter polar- 
ization measurements with photon energy, possibly providing a 
diagnostic for the maximum electron cutoff energy. 

2. Bremsstrahlung interactions and the creation of 
the photon distribution 

2.1. Stokes parameters and defining photon polarization 

The polarization state of incoherent radiation can be com- 
pletely described u sing four Stokes parameters dStokes|[l852t 
IChandraseldiarl[l960l) . The Stokes vector consists of these four 
parameters and takes the form of S = [I, Q,U,V]. The first 
Stokes parameter / is the normalised total intensity of the pho- 
ton beam, while Q/I, U/I and VI I will have values between -1 
and 1. The second and third normalised Stokes parameters are 
used to define linear polarization with 1 or -1 indicating that 
the beam/photon packet is completely polarized with the sign 
providing the direction of polarization. The fourth parameter is 
used to describe circular polarization but bremsstrahlung emis- 
sion in the solar corona, or Compton scattering only produces 
radiation that is linearly polarized. In order to produce circu- 
larly polarised radiation via bremsstrahlung, the spins of the ra- 
diating electrons need to be aligned and the magnetic field in 
the corona/chromosphere is obviously not strong enough for this 
alignment. This means that only the first three Stokes parame- 
ters are required and the fourth can be set to zero throughout our 



simulations. Generally in X-ray and gamma ray astronomy the 
polarization of radiation is measured using the degree of polar- 
ization (DOP) and the polarization angle which is the pre- 
ferred direction of the electric field. These are defined using the 
Stokes parameters as, 



DOP = 



Ve 2 + u 2 



and 



1 l-U 
¥ = - arctan 

2 l-Q 



(1) 



(2) 



where the angle *F is chosen to lie within the quadrant be- 
tween [-180°, 180°], so that arctan(±§) = +0, arctan = 

+ 180°, arctan (^) = -0 and arctan (fg) = -180°. The nega- 
tives introduced into equation (|2} ensure that a negative Q gives 
and a positive Q gives 90°. Hence with this definition, when 
*P = 0°, the observed radiation is polarized parallel to the radial 
direction at the solar disk and when *P = 90°, the radiation is 
polarized perpendicular to the radial direction. The opposite def- 
inition (i.e. 1/ = 90° - polarized parallel to the radial direction) 
would have been equally valid as long as definitions throughout 
the simulations are consistent. 

We notice that the Stokes parameters are also frame depen- 
dent and hence have to be upda ted by the use of rotation matrices 
dHovenier & van der Med 19831) when moving between different 
coordinate frames. The rotations used in our simulations when 
moving between the source frame and the scattering frame and 
vice versa are shown in appendices lAl and IB1 The DOP remains 
unchanged by the rotation but the polarization angle *P is mea- 
sured with respect to the new frame. 

For simplicity, the flare loop and hence the dominant direc- 
tion of electrons is always ass umed to lie parallel with the local 
solar vertical, i.e. no tilt (see lEmslie et al.l I2008L for non-zero 
tilt). This means that for the bremsstrahlung emission, the po- 
larization direction *P across the entire source (source assumed 
to be small in comparison with the solar disk) and at spatially 
resolved positions across the source can only ever equal 0° or 
90°, sin ce U from the brems strahlung emission is always close 
to zero (IBai & Ramatvl lT978). Since in the solar disk frame (the 
observed frame of the HXR source), Compton scattering can 
produce values of U other than 0, this means that the Compton 
scattered *P can have values other than 0° or 90°. This is true 
for the spatially resolved polarization of the albedo (and hence 
the total observed) source, while the spatially integrated albedo 
Stokes parameters again sum to produce *P values of either 0° or 
90°, for a non-tilted loop. 

2.2. Primary photon distribution 

The intensity of the photon distribution I(e,8) produced by 
bremsstrahlung for a chosen electron distribution F(E,/3) is 
given by 



I(e, 6) oc 



r*co r*2n r*7T 

I I I F(E,/3)o-,(E,e,Q)sm/3d/3d(bdE 

JE=e J<S>=0 J/3=0 



(3) 



where e is the photon energy, E is the electron energy, F(E,f3) 
is the target electron flux density differential in energy and 
crj(E, e, 0) is the total (averaged over all pol arization states) 
angle-dependent bremsstrahlung cross-section (lElwert & Haug 
119701: IBai & Ramatvl[T978T: iMassone et al.ll2004 . e [0°, 180°] 
is the photon polar emission angle measured from the local solar 



2 



Jeffrey and Kontar: Spatially resolved hard X-ray polarization in solar flares: effects of Compton scattering and bremsstrahlung 



vertical with = 0° directed away from the Sun. j3 e [0°, 180°] 
is the pitch angle of the emitting electrons velocity to the lo- 
cal magnetic field with ft = 0° also directed away from the Sun. 
<$> e [0°, 360°] is the corresponding electron azimuthal angle and 
&(J3, <t>, 9) is the angle between the plane of emission (at angle f3) 
and the plane of observation (at angle 9). The photon emission 
is described by ju = cos 9, where fi from to 1 corresponds to 
emission away from the Sun, and // from -0 to -1 corresponds to 
emission towards the solar surface. The photon emission angle 
fj. = cos 9 is related to the electron pitch angle ft by: 



cos = cos 9 cos /3 + sin 9 sin /3 cos <J). 



(4) 



Viewing the outward emission from fi = 0.0 -1.0 corresponds to 
observing the source at a selected heliocentric angle on the solar 
disk, i.e. fi = 0.0 would correspond to 90° and is equivalent 
to viewing a source sitting at the solar limb. The other Stokes 
parameters Q and U can be calculated in a similar manner: 



Q(e,9)°c I J I F(E,j3)crQ(E,e,®)sm/3d/3d®dE, 
JE=6 Jo=o J/3=0 

( 

E=e Jo=0 Jfl=0 



U(e, 9) oc 



F(E,p)a- u {E, e, 0) sin /3d/3d®dE, 



(6) 



with the only difference being the use of either <ti{E, e, 0), 
o- Q (E,e,@) or cru{E,e, 0). Q{e,9) and U(e,9) are then 
normalised between [-1,1] by dividing through by I(e,9). 
ctj(E, e, 0), ctq(E, e, 0) and cr u (E, e, 0) are the polariza- 
tion dependent cross-sec tions for bremsstrahlung taken from 
Gluckstern & Hull d!953h and also f ollowing the form used in 
Haujddl972lT and Emslie et al] d2008l) . are given by 



and 



07(E, e, 0) = cr ± (E, e, 0) + cr^E, e, 0), (7) 
<r Q {E, e, 0) = (crJE, e, 0) - cr^E, e, 0)) cos 20, (8) 



cru(E, e, 0) = (<t ± (E, e, 0) - cr n (E, e, 0)) sin 20, (9) 



where cr ± (E, e, 0) and <t\\(E, e, 0) are the perpendicular and par- 
allel components of the bremsstrahlung cross-section respec- 
tively. 

2.2.1. Electron distribution 

In our simulations the X-ray emitting electron distribution is 
chosen to have the form, 



F(E,p)ccE expf 



Av 2 



(10) 



The energy de pendence follows a power law as show n by ob- 
servations (see iHolman et al.ll201 UlKontar et al.ll201 ll as recent 
reviews) and the electron angular distribution is a Gaussian. The 
use of a Gaussian allows the angular anisotropy of the electron 
distribution to be easily controlled by a single parameter Av. The 
smaller the value of Av, the greater the proportion of the electron 
distribution, and hence the resulting photon emission directed to- 
wards the photosphere. 

In all simulations shown in this paper, we used a typical solar 
flare with 6 = 2 (spectral index of the mean electron spectrum 
and Av = 4.0 , Av = 0.5 or Av = 0.1 to describe various pitch 



angle distributions of electrons. Note that a target electron dis- 
tribution of 5 — 2 is produce d by an injected ele ctron distribu- 
tion of Sinj = 5 + 2 = 4- (see iHolman et al.ll2.Ql ll as a review). 
A Av = 4.0 electron distribution produces an approximately 
isotropic, unpolarised photon distribution, while the Av = 0.5 
and Av = 0. 1 electron distributions produce photon distributions 
with greater and greater beaming towards the photosphere. The 
resulting fluxes and polarizations of each component [primary 
towards the observer only (orange), albedo (blue) and the total 
observed emission (green)] plotted against emission/heliocentric 
angle yu = cos 9 are shown for each of the Av = 4.0, Av = 0.5 
and Av = 0.1 photon distributions in Figure Q] The resulting 
photon distributions in Figure Q] are shown for the energies of 
20-50 keV, where the albedo emission peaks. As expected, even 
though the Av = 0. 1 distribution has a greater downward beam- 
ing (and hence a smaller observer directed emission) than the 
Av = 0.5 distribution, within this low energy band the difference 
in anisotropy between the two primary distributions cannot be 
clearly seen, however differences between the two distributions 
can be clearly observed in the DOP with disk location over this 
energy range (second row Figure[T). 

A HXR source is modelled at a chosen height h in 
the chromosphere. For the HXR energy range, the source 
is about 1 Mm above the photosphere i n acco r dance with 
recent X-ray observations dKontar et al.1 120081: iPrato i 



20091: ISaint-Hilaire et all 120101: iMrozek & Kowalczukl 1201 



et al] 



iKontar et a l. 2010; Battaglia & Kontar 201 1). The source spatial 
extent is modelled using a two-dimensional circular Gaussian 
distribution, 



I(x, y) oc exp - 



2d 2 



y 

2d 2 



(11) 



where d is the size of the HXR source. Source sizes will be de- 
fined by the (Gaussi an) Ful l Width Half Maximum (FWHM), 
given by FWHM= V8 ln(2)d. 



3. Monte Carlo simulations of photon transport 



A Monte Carlo code developed bv lKontar & Jeffre"vl (2010) was 
modified to include polarization (see appendices A and B, for de- 
tails). The energy, angular and polarization properties of the in- 
put photon distributions for the MC code are determined via the 
chosen input electron distributions given by Equation [10] while 
the spatial properties are determined via Equation [TT] To input 
l(e,0), Q(e,9) and U(e,ff) (Equations [3] and |6]) into our MC 
simulations, distributions of energy e, angle 9 and polarization Q 
and U are numerically created from I(e, 9), Q(e, 9) and U(e, 9). 
Each photon azimuthal angle (f> is simply drawn from a uniform, 
random distribution between and 2n. To achieve the best statis- 
tics possible, we use 10° photons in every MC code run. Since 
the photon angular distribution is measured from the local solar 
vertical, photons with 9 < 90° escape from the solar atmosphere 
as the primary emission and photons with 9 > 90° will travel to- 
wards the photosphere, creating the albedo emission. Photons 
are allowed to move freely until they reach the photosphere, 
where they interact with free or bound electrons. In the simula- 
tions the photosphere is just defined by having a Hydrogen num- 
ber density of n H = 1.16xl0 17 cm' 3 dVernazza et al.ll98lh . Here 
the photons can no longer travel without interacting in some way 
with the photospheric material. Each photon moves a step-size 
before an interaction and this is calculated using ss = —l\n( step , 
where / is the photon mean free path and ( step e [0,1]. The 
photon mean free path is calculated by / = 1 /riH(r tota i, where 
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A;/ = 4.0 Af=0.5 Ai/=0.1 




fi=cosB /j,= cosd /i=cos5 

Fig. 1. Photon flux (top row) and DOP (bottom row) for upward primary (orange), albedo (blue) and total (green) components across 
the entire source for both the Av = 4.0 (1st column), the Av = 0.5 (2nd column) and the Av = 0.1 (last column) created photon 
distributions respectively, at the peak albedo energies of 20-50 keV, for different locations on the solar disk from the centre (pi = 1.0) 
to the limb (ji = 0.0). Note here that a negative DOP denotes that the direction of polarization *F across the source i s parallel to the 
radial direction. The DOP for Av = 4.0 photon distribution (near isotropic distribution) shows ~ the same result as iBai & Ra matv 
(Il978l) . as expected. 



Cwtai - o~c + o~a, the addition of th e Compton scat t ering cross- 
section (using the form shown in IBai & Ramatvl ([1978)) and 
photoelectric absorption cross section. This photon mean free 
path, /, is of the order 100 km. A photon interaction with the 
photospheric medium can either be by Compton scattering or 
by photoelectric absorption. For photons with energies less than 
~ 10 keV, photoelectric absorption is the more probable process 
while scattering dominates above ~ 10 keV (Figure [2]). When 
a photon is absorbed, it is removed from the simulations. For 
each photon, one of the two processes is chosen by calculat- 
ing the ratio of o>/cr to , fl /. Another random number f p! - c * is then 
sampled from a uniform distribution between and 1 . If the ra- 
tio is greater than then the photon is Compton scattered 
and if the ratio is less than (pick, then the photon is absorbed. 
The curvature of the Sun is included in these simulations and 
a photon exits the photosphere when it satisfies the condition 

z > zq - yjR 2 Q - x 2 -y 2 - Rq, where Rq is the radius of 

the Sun (taken to be 6.96 x 10 10 cm ~ 960"). The extent of 
the albedo patch is limited by properly modelling the curvature 
of the Sun. The photons are allowed to scatter multiple times 
until they exit the photosphere or are removed via absorption. 
Most photons will leave the photosphere during the first scatter, 
with subsequent scatterings producing less a nd less photons (th is 
is shown in terms of Green's functions in IKontar et al. (2006)). 
Photons that exit the photosphere with cos 6 > can then be 
collected into selected angular or energy bins corresponding to 
HXR sources sitting at any chosen heliocentric angle on/above 
the solar disk. In these simulations, source positions and sizes of 
each component (primary and albedo) and the total source are 
give n by the first and second moments of each intensity distribu- 
tion dKontar & Jeffrevl2010h . which are the mean (x, y) and vari- 
ance (var x , var,,) respectively. Though the primary distribution is 



initially Gaussian, the albedo (and hence total observed) distri- 
bution will have a complex shape that is no longer Gaussian. To 
quantify the sizes we use 'Gaussian' Full Width Half Maximum 
(FWHM) defined as 2.35 times the square root of the second- 
moment of the spatial distrib ution (not the actual F WHM of the 
complex distribution) (see IKontar & Jeffrey! [20101 for details). 
This allows a simple comparison with measurable parameters. 
All FWHM values in this paper are measured in this way. 



3. 1 . Compton scattering 

Above ~10 keV, Compton scattering is the dominant pro- 
cess in the photosphere. The polarization depe ndent differential 
Comp ton scattering cross-section is given bv iKlein & Nishina 
(1929), but the form used in these simulations is from lMcMastei 
d 196 II) and lBai & Ramatvl d!978l) 



dcr c 
dQ. 



1 



et) 



— + — 



•sin 2 6s 1 1 - Q cos 2<ps - U sin 20s 



))• 



(12) 

where ro = 2.82 x 10 13 cm is the classical electron radius, eo is 
the energy of the incoming photon, e is the energy of the outgo- 
ing photon, 6s is the polar scattering angle, <ps is the azimuthal 
scattering an gle and Q and U are the linear Stokes parameters 
respectively iMcMasterlll96lHBai & Ramatvlll978l) . The maxi- 
mum change in DOP occurs when 6s = 90° and no change in 
DOP occurs for a backscattering at 180°. The azimuthal scat- 
tering angle <ps also has a non uniform dependency on the in- 
coming polarization state (see Appendix IB1. In the simulations, 
the Klein-Nishina cross-section is multiplied by Z p i, olo = 1.18 to 
take account of elements higher than Hydrogen that are present 
within the photosphere. Z p i wto indicates the average atomic num- 
ber and the number of electrons per Hydrogen atom in the pho- 
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Fig. 2. Absorption cr a (black solid) and Compton cr c (red dashed) 
cross-sections multiplied by 10 24 e 3 (e in keV) at low energies 
below 10 keV. The absorption cross -section is calculated using 
photospheric element abundances by Asplun d et alj d2009t) . 



^observer 
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Fig. 3. Diagram of a single Compton scattering in the photo- 
sphere for three heliocentric angles of 0°, 45° and 90°. For a 
single scatter, the DOP tends to 100% as the scattering angle 
approaches 90°, producing a variation in polarization across the 
extent of the source. 



tosphere, and is given by 



-^photo 



£Z10 Az 

z 

10 A " 



(13) 



where Az is the log 10 abundance of an element with atomic num- 
ber Z relative t o Hydrogen while Ah — 12 is the log 1() abundance 
of Hydrogen dAsplund et al.l 12009). More details regarding the 
updating of photon parameters (6s ,£,<ps, Q, U) after Compton 
scattering are given in the appendices (A and B). 



3.2. Photoelectric absorption of X-ray photons 

For photons with energies below ~10 keV, photoelectric ab- 
sorption is the most probable photon interaction in the photo- 
sphere. The process of absorption is heavily dependent upon 
the composition of elements and the abundance of these ele- 
ments within the photosphere. Absorption was therefore mod- 
elled using the latest know n solar photospheric abundances taken 
from Asp lund et al.l d2009h . Absorption cross-section codes for 
the most important elements of H, He, C, N, O, Ne, Na, 
Mg, Al, Si, S, CI, At, Ca, Cr, Fe and Ni were adapted 
from Bal ucinska-Church & M cCammon (1992). Absorption is 
much less probable above 10 keV and for energies higher 
than 10 keV, the absorption cross-section was approximated 
by <T a (eo) °c e Q 3 . The probability of photoele ctric absorption 
is ass umed to be independent of polarization dPoutanen et al.l 
Il996l) . Only the angular distribution of the ejected electron is 
dependent upon the photon polarization, which is not mod- 
elled in our simulations. A comparison of the Compton scat- 
tering and absorption cross sections is shown in Figure 2 (with 
tr a and cr r multiplied by 10 24 e 3 (e in keV) for comparison 



with lMorrison & McCammorjjj 1983 )). Any differences between 
Figure [2] and Morrison & McCammo n! (Il983|) a re due to the 
newer element abundances (lAsplund et al.l 120091) and updated 
absorption cross-section codes (in particular Helium) being used 
in our simulations. 



4. Spatial distribution of X-ray polarization 

4.1. Single Compton scatter for an isotropic unpolarised 
source 

In order to demonstrate the spatial variation in polarization due 
to Compton scattering, the easiest example to consider is the 
albedo patch created by an initially isotropic, unpolarised point 
source at a height h above the photosphere (see Kontar & Jeffrey 
(2010)). For this example, the variation in po larization acros s the 
source can be described analytically by (McMaster 1961) and 
assuming no energy losses, 



DOP = 



1 - cos 2 6 S 
1 + cos 2 6 S 



where 



cos 6s = cos 6 cos(7r - 6i) + sin 6 sin(7r - #,) cos ( 



(14) 



(15) 



6s is the scattering angle, 6, is the emission angle measured from 
the local solar vertical, 6 is the heliocentric angle on the solar 
disk and <p is the angle measured in the solar disk plane. In equa- 
tion ( TBI , the scattering angle 6s determines the DOP. 6s is re- 
lated to distance r from the centre of the albedo patch by equa- 
tion (TT~5T > though r — h tan(7r - #,), and hence the DOP at any 
point across the albedo patch sitting at any heliocentric angle on 
the solar disk can be easily calculated. Note that for this simple 
example, equation <TT~4T > assumes that the energy difference be- 
tween the incoming and scattered photon is negligible, though 
this is only the case for low HXR energies of ~ 10 keV. The 
DOP for all photon energies is calculated using the T scatter- 
ing matrix, which is used in t he simulations (equation dB.3t in 
appendix|B]i (McM asterll 19611) . 

Figure [3] shows a cartoon of this single scattering from an 
isotropic point source for three different heliocentric angles of 
0°, 45° and 90°. For a source located exactly above the solar 
centre (0°) at a height h, the resulting variation across the pho- 
tospheric albedo patch is radially symmetric. As the radial dis- 
tance r from the source centre increases, the scattering angle 6s 
of any observed radiation will decrease from 180° towards 90°, 
causing the DOP to grow from 0% to ~ 100%. Radiation scat- 
tered in the photosphere at a location directly below the source 
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(180° backscatter) and emitted towards the observer will experi- 
ence no change in its DOR This statement is true for sources at 
any heliocentric angle 0, but projection effects at angles 0, > 0° 
will create an asymmetry in the polarization pattern along the 
radial direction, whereas the polarization pattern in the perpen- 
dicular to radial direction always remains symmetrical. The de- 
scribed pattern can be seen in Figure [4] which shows the polar- 
ization maps for a single Compton scattering at four disk loca- 
tions ranging from the solar centre to the limb, along Y = 0" 
atX = 214", 543", 750", 936". These locations are equivalent 
to fi = 0.97,0.82,0.62,0.22 (0, = 14°, 35°, 52°, 77°) and de- 
note the approximate positions of the total source, which are 
shift ed from the primar y position due to the albedo component 
dKontar & Jeffrey! 12010). The polarization across the source at 
any location on the solar disk can always be measured with re- 
spect to the radial connecting the solar disk centre and the cen- 
tre of the source. Therefore due to the symmetry of the prob- 
lem, source locations at Y = 0" considered in this paper can 
straightforwardly be applied to any solar disk location. In Figure 
|H the dotted blue ellipse denotes the FWHM of the albedo com- 
ponent, the solid green ellipse denotes the FWHM of the total 
observed source and the orange, blue and green asterisks indi- 
cate the centroid positions of primary, albedo and total observed 
sources respectively. The polarization angle *P follows a steady 
pattern across all sources. *F is always at an angle tangential to 
the line connecting the desired position and the location of the 
source centroid position. Hence at disk centre locations the *P 
pattern is also symmetrical across the source. 

Figure [5] shows the same albedo polarization maps as in 
Figure |4]but for multiple Compton scatterings. The overall pat- 
tern for the DOP and *P are preserved but due to multiple scatter- 
ings the overall DOP across all locations over the albedo patch 
has decreased. A single scattering DOP of ~ 100% near the edge 
of the source has been reduced to ~ 50% by multiple scatterings. 
All other simulations shown below are for multiple Compton 
scatterings. 

The albedo pattern from a primary source at a greater height 
than 1 Mm (say from a coronal source) should produce the 
same albedo polarization pattern but over a much greater area 
in the photosphere. Plotting the polarization pattern at a higher 
or lower energy than the peak albedo range of 20-50 keV would 
obscure the albedo polarization pattern due to lowering albedo 
fluxes at these energies. 



4.2. Anisotropic source at a height ofh=l Mm (1.4") and 
size of 5" 

For a chosen chromospheric HXR source height of h — 1 
Mm (1.4") and a primary source size of FWHM~ 5", simula- 
tions were performed for all three photon distributions created 
by the Av = 4.0, Av = 0.5 and Av = 0.1 electron distribu- 
tions. This height was chosen to match chromos p heric H XR 
sourc e height measurements (pCpntar et ail 120081 : iPrato et aT 



200 91 ISaint-Hilaire et all 1201 Ot lMrozek& Kowalczukl |201( 



Kontar etal J 120 lOt iBattaglia & Kontarll201lh . All the results 
shown here are for the energy range of 20-50 keV, where albedo 
emission peaks, producing the largest distortion to the primary 
component but the best opportunity for the detection of the 
albedo component. 

Figures [6] [9] and [12] each plot the resulting polarisation 
maps for four HXR sources (resulting from the primary and 
albedo components) created by the Av = 4.0, Av = 0.5 and 
Av = 0. 1 electron distributions respectively. As with Figures [4] 
and [5] each figure plots four HXR sources positioned at fi ~ 



0. 97.0.82.0.62.0.22. In Figures [6] [9] and [[2] the dotted ellipses 
denote the FWHM of the total source (green), the primary source 
(orange) and the albedo source (blue) and the correspondingly 
coloured asterisks denote the (x,y) centroid position of the to- 
tal source and the primary and albedo components respectively. 
Figures [7] [10] and [13] plot intensity, I (top row), DOP (middle 
row) and *P (bottom row) along the radial direction X centred 
at Y = 0" (across a of bin width= 2") for each of the maps in 
Figures[6l[9]and[T2]. Figures[8][II]and[T3plot the DOP (top row) 
and *P (bottom row) along the perpendicular to radial direction 
Y centred at X= 213.6", 543.9", 750.1", 936.6" (again across a 
bin width= 2" ) for each of the maps in Fi gures l6l l9l and [T2l 

4.3. Quasi-isotropic distribution with Av = 4.0 

The photon distribution produced by the Av = 4.0 electron dis- 
tribution is approximately unpolarised and isotropic. Therefore, 
both spatially integrated and spatially resolved polarization mea- 
surements across the primary source at all locations on the solar 
disk produce DOP- 0% and ¥ = 0° (radial) at 20 - 50 keV. 
The albedo component produces asymmetrical DOP and 4* vari- 
ations along the source radial direction (Figure [7]) while along 
the source perpendicular to radial direction (Figure[8]), variations 
in albedo DOP and *P are approximately symmetrical, since the 
centroid positions of the primary and albedo components always 
coincide in the perpendicular to radial direction. The simulated 
HXR sources plotted in Figure [6] have a finite source size of 
~ 5". Compared with a point source (Figure |5), this produces 
two main differences: i) photons leave the source from differ- 
ent positions above the photosphere and ii) for certain distri- 
butions and disk locations, the primary polarization will dom- 
inate over the extent of the primary source. The first property re- 
duces the DOP at all points across the source, compared with the 
albedo patch created by a point source. Due to the second prop- 
erty, the polarization variation caused by albedo may be slightly 
masked by the primary component within the source FWHM, 
especially for cases where the primary component is dominant, 

1. e. isotropic or near isotropic distributions. 

For the quasi-isotropic Av = 4.0 distribution, the primary 
component is the dominant component at all four disk locations, 
with the albedo contribution falling as the source location moves 
towards the limb (Figure |7J- first row). Hence, the primary com- 
ponent dominates within the FWHM of the total source, while 
the albedo component dominates after this boundary. The second 
and third rows of Figure [7] demonstrate common radial trends 
in DOP and T*, not only across the extent of each individual 
source but also between sources at different disk locations. For 
a quasi-isotropic distribution, at a particular disk location (other 
than the disk centre), the highest DOP along the radial direc- 
tion is observed at the disk-centre-side of the source (where the 
albedo dominates). This falls to approximately zero within the 
FWHM of the total source (where the primary dominates) and 
then increases again towards the limb-side of the source (where 
the albedo again dominates), but always remains lower than the 
DOP at the disk-centre-side. Comparing the four disk locations, 
the DOP at all points along the radial direction of a single source 
decreases as the source location nears the limb. In the radial di- 
rection, spatially resolved DOP can achieve values as high as 
~ 30% at disk centre locations. The albedo component produces 
the distinctive T* variation shown in Figure [6] Along the radial 
direction, *¥ = 90° at the disk-centre-side of the source, falls to 
zero within the FWHM extent and then rises again at the limb- 
side of the source. Along the perpendicular to radial direction of 
a single source (Figure [8), DOP (first row) and the magnitude of 
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Fig. 4. Albedo polarization maps for an isotropic, unpolarised point source sitting above the photosphere at four different radial 
locations of X = 214", 543", 750", 936" at Y = 0" (corresponding /j = 0.97,0.82,0.62,0.22) after a single Compton scatter in the 
photosphere. All results are shown at the peak albedo energies of 20-50 keV. The length of each red arrow indicates the DOP and 
the direction of each arrow depicts the polarization angle ¥ within the chosen plotting bin. The solar radial direction (or X axis for 
this case) is defined as the *F = 0° position. An arrow length of 2" corresponds to a maximum DOP of 100%. The green and blue 
ellipses give the FWHM of the total and albedo sources respectively, while the green, blue and orange asterisks give the centroid 
position of the total, albedo and primary sources. 
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Fig. 5. Albedo polarization maps as in Figure [4] but for the case of multiple Compton scatterings in the photosphere. Multiple 
scatterings have 

acted to decrease the DOP at all points across each source. 



*P (second row) are symmetrical due to no projection effects and 
the centres of the primary and albedo components always coin- 
ciding. As with the radial direction, in the perpendicular to radial 
direction, spatially resolved DOP can achieve values as high as 
~ 30%. While the magnitude of VP across the source is symmetri- 
cal at each disk location in the perpendicular to radial direction, 
*P itself behaves as an odd function along Y, with a 180° ro- 
tational symmetry about the source centre, increasing from the 
radial at the upper source edge to IT* > 0°| and then back to the 
radial direction at the lower source edge. 

4.4. Beamed electron distributions Av = 0.5 and Av = 0.1 

All plots for the photon distribution created by the Av = 0.5 
electron distribution are shown in Figures (l9lfTT1l. while Figures 
( fT2lfT4l i show all plots for the photon distribution created by the 
Av = 0.1 electron distribution. Comparison of Figures [6] [9] and 
[T2]demonstrates that increased beaming towards the photosphere 
produces smaller, more concentrated and intense albedo patches. 

For the Av = 0.5 distribution, Figure [10] plots the radial in- 
tensity, I (top) (along X at Y — 0"). The first two disk locations 
are albedo dominated, the third disk location has approximately 
equal contributions from the primary and albedo components 



and only the disk location closest to the limb is primary dom- 
inated. The primary DOP can rise as high as ~ 20% at the limb 
and the primary \P is radial at all locations. Figure[13](top) plots 
radial intensity slices (along X at Y = 0") for the Av = 0.1 dis- 
tribution. As expected, the first three disk locations are albedo 
dominated and the primary DOP can reach ~ 40% at the limb. 
Again, the primary \P is radial at all disk locations. 

As with the quasi-isotropic Av = 4.0 distribution, common 
trends can be observed across individual sources at particular 
disk locations and between disk locations for both the Av = 0.5 
and Av = 0.1 distributions. More importantly for observations 
and anisotropy deduction purposes, trends between each of the 
three simulated distributions (Av = 4.0, Av = 0.5 and Av = 0.1) 
can be observed, along the radial and perpendicular to radial di- 
rections, at any chosen disk location. The most notable trends 
are observed in the radial (X) direction, and it is these trends that 
may help deduce something about the anisotropy of the photon 
distribution for a HXR source sitting at a given disk location. 
Trends can be observed at all disk locations, but in this exam- 
ple the patterns are most noticeable in the second and third disk 
locations plotted. In both of these locations, the disk-centre-side 
DOP falls with increased beaming while the limb-side DOP rises 
with increased beaming. 
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Fig. 6. Total X-ray brightness and polarization (arrows) maps for the photon distribution created by the near isotropic Av = 4.0 elec- 
tron distribution for a 5" primary source. The resulting total sources sit at four disk locations of X = 213.6", 542.9", 750.1", 936.6" 
at Y — 0" (corresponding to ft ~ 0.97, 0.82, 0.62, 0.22). The green, blue and orange ellipses and asterisks give the FWHM and 
centroid positions of the total, albedo and primary sources respectively. The primary component is dominant at all disk locations 
and hence the albedo polarization pattern at small radii (over the extent of the 5" primary source) is obscured slightly by the fiat, 
constant polarization of the primary source. 



Comparing the third disk location (for example) in Figures 
171 [TOl and 1 131 shows how the radial DOP at the limb-side of the 
source rises with increased beaming from ~ 2% for the Ay = 4.0 
distribution to ~ 18% for the Av = 0.5 distribution to ~ 30% for 
the Av = 0.1 distribution. The radial DOP at the disk-centre- 
side of the source falls with increased beaming from ~ 18% 
(Av = 4.0) to ~ 14% (Av = 0.5) to ~ 4% (Av = 0.1). The 
polarization angle *P also produces similar patterns with chang- 
ing anisotropy. A clear example of this can be observed by com- 
paring the second disk location plotted in Fi gure s 171 [TOl and [T3l 
Along the radial direction, disk-centre-side *P generally stays 
at SF = 90° for all photon anisotropies, while the outer limb- 
side *P falls significantly with increased beaming, from *P = 60° 
(Av = 4.0) to ¥ = 20° (Av = 0.5) to Y = 0° (Av = 0.1). 

Therefore, the DOP and SP patterns are a clear indication of 
how spatially resolved polarization measurements could be used 
to determine the beaming of the photon distribution. It should be 
noted that a (near) disk-centre source produces a slightly differ- 
ent trend in radial DOP with increasing photon anisotropy. The 
DOP at the disk-centre-side remains approximately the same for 
all photon anisotropies while the DOP at the limb-side falls with 
increased beaming (this is the opposite trend to other disk loca- 
tions). 

Comparing each disk location along the perpendicular to ra- 
dial direction (Y) in Figures [8] QT| and [14] shows that greater 
beaming increases the DOP over the whole extent of the source 
at any given location (except at the disk centre where the spa- 
tially resolved polarization along Y is approximately the same 
for all three distributions). This spatial increase is most notice- 
able at limb locations where from the source centre to the source 
edge, DOP increases from ~ 0% to ~ 20% (Av = 4.0), from 
~ 20% to ~ 40% (Av = 0.5) and from ~ 40% to ~ 55% 
(Av = 0.1). 

4.5. High energy cutoff in the electron distribution and HXR 
polarization 

Spatially integrated polarization is dependent on the highest 
energy in the electron distribution (the high cutoff energy) 
dHeristchill987l) . When calculating spatially integrated polariza- 
tion, equation ([TJ reduces to DOP - Q and equation (f2]i reduces 
to *P = 0.5 arctan(^) as U sums to zero for a single measure- 



ment across the entire source. A negative DOP indicates that the 
polarization angle is parallel to the radial direction (*F = 0°), 
while a positive DOP indicates that the polarization angle is per- 
pendicular to the radial direction Q¥ = 90°). 

For the three electron distributions of Av = 4.0, Av = 0.5 
and Av = 0.1, simulations were run with two high cutoff elec- 
tron energies of E cut0 ff = 500 keV and E cut0 ff = 2 MeV. Figures 
Q3][T6]and[T7]plot the spatially integrated polarization across the 
total source (green) and the primary source only (orange) against 
photon energy e at four disk locations pt € [0.20 - 0.25], [0.60 - 
0.65], [0.80 - 0.85], [0.95 - 1 .00] for each of the above distribu- 
tions respectively. The important property to observe here is not 
the magnitude but the sign of the DOP (or whether ¥ = 0° or 
¥ = 90°). 

Using an electron distribution with Av = 4.0 (FigurefTBIl and 
a cutoff energy of E cut0 ff = 500 keV produces a photon dis- 
tribution with a negative DOP at all photon energies and disk 
locations, while using the same distribution with a cutoff energy 
of E cutoff = 2 MeV creates a photon distribution where the DOP 
changes from negative to positive at ~ 100 - 200 keV at all disk 
locations. During bremsstahlung, in order to conserve energy, 
electrons with higher energies will scatter through larger angles. 
When a photon is scattered through a large angle, its polariza- 
tion is more likely to be directed perpendicular to the plane of 
emission (*P = 90°) rather than parallel to the plane of emission 
0? = 0°). Therefore, a change in the direction of polarization 
(from *F = 0° to *P = 90°) indicates the presence of higher ener- 
gies in the electron distribution, greater than ~ 1 MeV. 

As the beaming of the electron distribution increases (the 
Av = 0.5 and Av = 0.1 distributions), we see that the above state- 
ment does not hold and it becomes more likely that both electron 
distributions with cutoff energies of either E cuto ff = 500 keV 
and E cuto ff = 2 MeV will produce photons with *P = 90°. For 
the Av = 0.5 distribution (Figure [16), we see that as the source 
moves towards the solar centre, the photon distribution created 
by the E CMO ff = 500 keV also produces photons at higher ener- 
gies (again ~ 100-200 keV) with *P = 90°. For the very beamed 
Av = 0.1 distribution (Figure [TTJ, both the E CUW ff = 500 keV 
and Ecutoff - 2 MeV distributions produce high (~ 100 - 200 
keV) energy photons with ¥ = 90° at all disk locations. 

Therefore using the direction of *P as an indicator for high 
energies in the electron distribution becomes less and less useful 
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Fig. 7. Radial slices through each of the sources in Figure|6]for the Av = 4.0 created photon distribution. Each of the radial slices are 
taken along X at Y = 0" for the intensity, /, DOP and polarization angle *P at the four source locations. As before, orange=primary, 
blue=albedo and green=total. The green dash-dot lines denote the centroid position and the FWHM of the total observed source 
while the orange and blue lines denote the centroid positions of the primary and albedo sources respectively. Slices along the radial 
direction for the Av = 4.0 distribution show a clear DOP and *F pattern at all four disk locations plotted. 



as the beaming of the photon (electron) distribution increases. 
The increased beaming causes lower and lower energy photons 
to scatter at larger angles, especially at locations closer to the 
solar centre, hence producing photons with ¥ = 90° in the 
Ecutoff = 500 keV distributions. Note that this method is only 
useful when the beaming of the photon distribution has already 
been determined. 



5. Discussion and conclusions 

The simulation results show that Compton backscattering pro- 
duces a clear and distinct albedo polarization pattern across an 
HXR source at energies of 20-50 keV. Trends can be observed 
for both of the measured polarization parameters, DOP and *P, 
and are clear indications of the existence of an albedo component 
in comparison with the constant, radial polarization of the pri- 
mary emission. This means that spatially resolved polarization 
can be used to probe structure within HXR footpoint sources, 
helping us distinguish between the bremsstrahlung source and 
the albedo source. More importantly, at a single disk loca- 
tion, spatially resolved DOP and *P exhibit clear variations with 
changing photon directivity, along both the radial and perpen- 
dicular to radial directions and can be used to determine the 



anisotropy of the electron distribution. Therefore, to take advan- 
tage of these properties requires reliable polarization measure- 
ments with an angular resolution of ~ 5" - 10" over the peak 
albedo energies of ~ 20 - 50 keV. 

The simulations also suggest that for ~ isotropic HXR 
sources, spatially integrated polarization angle measurements, 
from low to high energies, with consideration of the changes 
due to albedo, could help indirectly infer the highest energy 
in the electron distribution. For near isotropic sources im- 
plied by RHESSI X-r ay observations (Kaspa rova et alj [2007; 
iKontar & Brownll2006l) . changes in ¥ from the radial Q¥ = 0°) 
to the perpendicular to radial direction Q¥ = 90°) may help indi- 
cate the presence of high energy electrons (~> 1 MeV) present 
in the electron distribution, if the photon anisotropy is known. 
Changes in spatially integrated DOP measurements, from low 
to high energies, will also help determine the anisotropy of the 
photon distribution. 

Currently, when observing solar flares we do not have the 
instrumentation required to produce spatially resolved polariza- 
tion measurements. It is doubtful whether near future missions 
such as t he Gamma - Ray I mager/Polarimeter for Solar Flares 
(GRIPS) dShih et all [20091) may have the capability to mea- 
sure polarization over albedo energies, even though it should be 
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Fig. 8. Perpendicular to radial slices through each of the sources shown in Figure [6] for the Av = 4.0 distribution. Each of the 
perpendicular to radial slices are taken along Y at X = 213.6", 542.9", 750.1", 936.6" for the DOP and polarization angle The 
lines and colours are as in Figures [6] The DOP and the magnitude of the polarization angle *F remain symmetrical along the 
source perpendicular to radial direction. 



20-50 keV 




205 210 215 220 535 540 545 550 740 745 750 755 930 935 940 945 

X [orcsec] X [arcsec] X [arcsec] X [arcsec] 

Fig. 9. Total X-ray brightness and polarization maps for the photon distribution created by the Av = 0.5 electron distribution 
for a 5" primary source. The resulting total sources sit at four disk locations of X — 213.1", 542.6", 750.0", 936.5" at Y = 0" 
(corresponding to /j ~ 0.97, 0.82, 0.62, 0.22). The green, blue and orange ellipses give the FWHM of the total, albedo and primary 
sources respectively. Green, blue and orange asterisks give the centroid position of the total, albedo and primary sources. For disk 
locations closer to the solar centre, the albedo component is the dominant component over the primary component. 



able to measure polarization across 12.5". The best polarization 
measurements are likely to be over the rang e of 150-650 keV 
with a ~4% minimum detectable polarization dShih et al.l l2009). 
At energies greater than ~100 keV, the albedo flux drops off 
steeply, so it is very unlikely that any albedo component could 
be detected at these energies and the polarization across the 
observed HXR source would only be from the bremsstrahlung 
emission. Therefore, good polarization measurements at 150- 
650 keV from flares with high fluxes will give us a direct mea- 
surement of the primary component and DOP/T measurements 
at these energies may be used to infer the high energy cutoff in 
the electron distribution, for relatively isotropic distributions. 
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Appendix A: Compton scattering and new 
scattering angles and energies 

In the Monte Carlo simulations, when a Compton scattering oc- 
curs, the properties of the outgoing photon: energy e, polar scat- 
tering angle 6s , azimuthal scattering angle (f>s and the two linear 
polarisation parameters Q and U need to be updated. New polar 
scattering angles 6s for each photon can be found by integrating 
the Klein-Nishina differential cross section over <ps to produce a 
differential cross-sect ion that is only depende nt on e and 6s . This 
is the Klein-Nishina dKlein & Nishinall 19291) differential cross- 
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Fig. 10. Radial slices (along X) through Y = 0" for the intensity, /, the DOP and ¥ for each of the sources in Figure [9] for the 
Av = 0.5 distribution. Again, orange=primary, blue=albedo and green=total. The green dash-dot lines denote the centroid position 
and the FWHM of the total observed source while the orange and blue lines denote the centroid positions of the primary and albedo 
sources respectively. 



section for a completely unpolarised isotropic beam and is given 
by, 



do- c (0s,e) = 1 
dQ 2' 

















f-5- 











(A.l) 



Since MC simulations operate by drawing numbers randomly 
from a given distribution, this means that 8$ can be easily found 
by matching each value of 6s e [0°, 180°] with a number gg e 
[0, 1] for every value of e using, 



(t 



_2n r 
v c Jo 



dcr c {6 s , e) 
dQ 



d0 s , 



(A.2) 



where cr c is the total Compton scattering cross-section. New val- 
ues of 6$ are simply drawn at each scattering using the photon 
energy before scattering and a random, uniform number between 
Oandl. 

Once the new scattering angles &s are obtained then the new 
photon energy e can be easily found using, 
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(A.3) 



Appendix B: Compton scattering and photon 
polarization states 

If the photon distribution is completely isotropic and unpolarised 
then the azimuthal scattering angle <p$ can just be sampled from 



a uniform distribution between and 2n, but this is not true for 
the more general polarization dependent case. The probability 
density function of obtaining a value of (ps between <ps and <p$ + 
dips can be described by: 

1 do- c (e,6 s ,<t>s)/dCl 



P(<Ps) = 



In do- c (e,6 s )/dQ. 



1 f + ^ - sin 2 6» s (l - Qcos2(f> s - t/sin20 s ) 



2n 



a + - sin 2 6 S 



with the maximum value of this function given by: 

I f + j--sm 2 e s (l- Ve^Tc^) 



Ms) = 



2n 



(2 + £- - sin 1 



!) 



(B.l) 



(B.2) 



Firstly, a value of 4>s is sampled between and 27r. The condition 
that P(<ps) < P maxims) is men use d t° accept a value of <ps and 
provides a method for sampling values of <ps for each photon, 
using the new values 6s and e already calculated for each pho- 
ton/photon packet. This method is repeated until the condition 
is satisfied for each photon and each photon is provided with an 
azimuthal scattering angle (ISalvat et al.ll2008l) . 

Due to Compton scattering, the Stokes parameters have 
to be updated using the scattering matrix T (lMcMasterlll%Tl 
iBai & Ramat" " 



T(e,6 s ) = 



( + a + sin 2 
sin 6, 



cos 



sin 

2, 





6 S + 1 
2cos# t 



(B.3) 
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Fig. 11. Perpendicular to radial slices through each of the sources shown in Figure [9] for the Av = 0.5 distribution. Each of the 
perpendicular to radial slices are taken along Y at X = 213.1", 542.6", 750.0", 936.5" for the DOP and polarization angle W. The 
lines and colours are as in Figures [9] [10] 




Fig. 12. Total X-ray brightness and polarization maps for the photon distribution created by the Av = 0.1 electron distribution 
for a 5" primary source. The resulting total sources sit at four disk locations of X — 213.0", 542.5", 749.9", 936.4" at Y = 0" 
(corresponding to /j ~ 0.97, 0.82, 0.62, 0.22). The green, blue and orange ellipses give the FWHM of the total, albedo and primary 
sources respectively. Green, blue and orange asterisks give the centroid position of the total, albedo and primary sources. At all disk 
locations the albedo component is the dominant component over the primary component. 



Before a scattering, the Stokes parameters have to be rotated by 
<ps so that they are defined relative to the plane of scattering. The 
Stokes parameters before scattering are rotated by the rotation 
matrix M(<p s ) given by: 

no o ^ 

M((/> s )= cos 20s sin 2^5 (B.4) 
,0 - sin 2<ps cos 2<p$ , 

After a scattering, the Stokes parameters have to be rotated again 
so that they are defined relative to the starting position of the 
source. The Stokes parameters after scattering are rotated by the 
rotation matrix M(— T). 



where 



M(-F) 



cos r = +- 



10 o ^ 

cos -2r sin -2F 
- sin -2F cos —2T , 

w' — cos G s w 

Vl - COS 2 I s Vl - W 2 ' 



(B.5) 



(B.6) 



w and w' are the current and previous z direction cosines respec- 
tively dHovenier & van de~ Mee1[l983h . The ± in equation ( IB. 61 ) 
is present due to the negative sign being used when n < <p s <2n 
and the positive sign for < <f>s < n. The angle F is shown 



in Figure IB. II T is the angle between the scattering plane and 
the normal plane in the frame of the source. Therefore during 
a Compton scattering the order of the rotations on the Stokes 
vector [IQU] is M(^ s )T{e,6 s )M{-Y). 
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Fig. 13. Radial slices (along X) through Y = 0" for the intensity, /, the DOP and *T for each of the sources in Figure [12] for the 
Av = 0.1 distribution. Again, orange ^primary, blue=albedo and green=total. The green dash-dot lines denote the centroid position 
and the FWHM of the total observed source while the orange and blue lines denote the centroid positions of the primary and albedo 
sources respectively. 
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Fig. B.l. The position of the photon before scattering (blue) and 
after scattering (green) and the angle T that determines the fi- 
nal rotation of the Stokes parameters back into the frame of the 
source from the scattering frame. 



Emslie, A. G., Bradsher, H. L., & McConnell, M. L. 2008, ApJ, 674, 570 
Gluckstern, R. L. & Hull, M. H. 1953, Physical Review, 90, 1030 
Haug, E. 1972, Sol. Phys., 25, 425 
Henoux, J. C. 1975, Sol. Phys., 42, 219 
Heristchi, D. 1987, ApJ, 323, 391 

Holman, G. D., Aschwanden, M. J., Aurass, H., et al. 201 1, Space Sci. Rev., 159, 
107 

Hovenier, J. W. & van der Mee, C. V. M. 1983, A&A, 128, 1 
Kasparova, J., Kontar, E. P., & Brown, J. C. 2007, A&A, 466, 705 
Klein, O. & Nishina, T. 1929, Zeitschrift fur Physik, 52, 853 
Kontar, E. P. & Brown, J. C. 2006, ApJ, 653, L149 

Kontar, E. P., Brown, J. C, Emslie, A. G., et al. 201 1, Space Sci. Rev, 159, 301 
Kontar, E. P., Hannah, I. G., Jeffrey, N. L. S., & Battaglia, M. 2010, ApJ, 717, 
250 

Kontar, E. P., Hannah, I. G., & MacKinnon, A. L. 2008, A&A, 489, L57 
Kontar, E. P. & Jeffrey, N. L. S. 2010, A&A, 513, L2+ 

Kontar, E. P., MacKinnon, A. L., Schwartz, R. A., & Brown, J. C. 2006, A&A, 
446, 1157 

Langer, S. H. & Petrosian, V. 1977, ApJ, 215, 666 
Leach, J. & Petrosian, V. 1983, ApJ, 269, 715 

Lemen, J. R., Chanan, G. A., Hughes, J. P., et al. 1982, Sol. Phys., 80, 333 
Massone, A. M., Emslie, A. G., Kontar, E. P., et al. 2004, ApJ, 613, 1233 
McConnell, M. L., Ryan, J. M., Smith, D. M., et al. 2007, in Bulletin of the 

American Astronomical Society, Vol. 38, American Astronomical Society 

Meeting Abstracts #210, 211-+ 
McConnell, M. L., Smith, D. M., Emslie, A. G., et al. 2004, Advances in Space 

Research, 34, 462 
McMaster, W. H. 1961, Reviews of Modern Physics, 33, 8 
Morrison, R. & McCammon, D. 1983, ApJ, 270, 119 

Mrozek, T. & Kowalczuk, J. 2010, Central European Astrophysical Bulletin, 34, 
57 

Muleri, E, Bellazzini, R., Brez, A., et al. 2009, Journal of Instrumentation, 4, 
11002 

Nakada, M. P., Neupert, W. M., & Thomas, R. J. 1974, Sol. Phys., 37, 429 



13 



Jeffrey and Kontar: Spatially resolved hard X-ray polarization in solar flares: effects of Compton scattering and bremsstrahlung 




Fig. 14. Perpendicular to radial slices through each of the sources shown in Figure [12] for the Av = 0.1 distribution. Each of the 
perpendicular to radial slices are taken along Y at X = 213.0", 542.5", 749.9", 936.4" for the DOP and polarization angle The 
lines and colours are as in Figures [12] [13] 
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Fig. 15. Changes in flux (top row) and spatially integrated polarization (bottom row) for the total source (green) and primary (orange) 
and albedo (blue) components with photon energy e for the Av = 4.0 created photon distribution with E cuto ff = 500 keV (solid 
lines) and E cuto ff = 2 MeV (dashed lines). Note that the four disk locations shown are the same as those used in Figure[6] Though 
there are small changes in the magnitude (DOP) of the polarization, at high photon energies (~ 200 keV), the polarization direction 
of the 2 MeV distribution changes from *P = 0° (negative) to *P = 90° (positive), while ¥ for the 500 keV distribution always stays 
negative. 
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Fig. 16. Changes in flux (top row) and spatially integrated polarization (bottom row) for the total source (green) and primary (orange) 
and albedo (blue) components with photon energy e for the Av = 0.5 created photon distribution with E cuto ff = 500 keV (solid 
lines) and E cuto ff = 2 MeV (dashed lines). 
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Fig. 17. Changes in flux (top row) and spatially integrated polarization (bottom row) for the total source (green) and primary (orange) 
and albedo (blue) components with photon energy e for the Av = 0.1 created photon distribution with E cuto ff = 500 keV (solid 
lines) and E cut0 ff = 2 MeV (dashed lines). 
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